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Abstract 

We derive a closed-form solution for the price of an average price as well as an average strike 
geometric Asian option, by making use of the path integral formulation. Our results are compared 
to a numerical Monte Carlo simulation. We also develop a pricing formula for an Asian option 
with a barrier on a control process, combining the method of images with a partitioning of the set 
of paths according to the average along the path. This formula is exact when the correlation is 
zero, and is approximate when the correlation increases. 
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I. INTRODUCTION 



Since the beginning of financial science, stock prices, option prices and other quantities 
have been described by stochastic and partial differential equations. Since the 1980s how- 
ever, the path integral approach, created in the context of quantum mechanics by Richard 
Feynman has been introduced to the field of finance j^yj]. Earlier, Norbert Wiener (4], 
in his studies on Brownian motion and the Langevin equation, used a type of functional 
integral that turns out to be a special case of the Feynman path integral (see also Mark 
Kac jlj], and for a general overview see Kleinert js] and Schulman Q). The power of path- 



integration for finance 



:iiJs- 



14|) lies in its ability to naturally account for payoffs that are 



path-dependent. This makes path integration the method of choice to treat one of the most 
challenging types of derivatives, the path-dependent options. Feynman and Kleinert [l5 | 
showed how quantum-mechanical partition functions can be approximated by an effective 
classical partition function, a technique which has been successfully applied to the pricing 
of path-dependent options (see ref. 6j and references therein, and Refs. 12|, ll6] for recent 
applications) . 

There exist many different types of path-dependent options. The two types which are 
considered in this paper are Asian and barrier options. Asian options are exotic path- 
dependent options for which the payoff depends on the average price of the underlying asset 



during the lifetime of the option 
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17- 



19) . One distinguishes between average price and 



average strike Asian options. The average price Asian option has been treated in the context 
of path integrals by Linetsky [2^]. The payoff of an average price is given by max(SV — K, 0) 
and max(i^ — St, 0) for a call and put option respectively. Here K is the strike price and 
St denotes the average price of the underlying asset at maturity T. St can either be the 
arithmetical or geometrical average of the asset price. Average price Asian options cost 
less than plain vanilla options. They are useful in protecting the owner from sudden short- 
lasting price changes in the market, for example due to order imbalances 21] ■ Average strike 
options are characterized by the following payoffs: max(S , r — St, 0) and max(Sr — St, 0) for 
a call and put option respectively, where St is the price of the underlying asset at maturity 
T. Barrier options are options with an extra boundary condition. If the asset price of 
such an option reaches the barrier during the lifetime of the option, the option becomes 
worthless, otherwise the option has the same payoff as the option on which the barrier has 
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been imposed, (for more information on exit-time problems see Ref. 22J and the references 
therein) 

In section [III we treat the geometrically averaged Asian option. In section III Al the asset 
price propagator for this standard Asian option is derived within the path integral framework 



in a similar fashion as in Ref. 



201 ] for the weighted Asian option. The underlying principle of 



this derivation is the effective classical partition function technique developed by Feynman 
and Kleinert [l^ . In section III Bl we present an alternative derivation of this propagator 
using a stochastic calculus approach. This propagator now allows us to price both the 
average price and average strike Asian option. For both types of options this results in a 
pricing formula which is of the same form as the Black-Scholes formula for the plain vanilla 



option. Our result for t 
found in the literature 



re option price of an average price Asian option confirms the result 
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23j . For the average strike option no formula of this simplicity 



exists as far as we know. Our derivation and analysis of this formula is presented in section 
III CI where our result is checked with a Monte Carlo simulation. In section II III we impose a 
boundary condition on the Asian option in the form of a barrier on a control process, and 
check whether the method used in section [III is still valid when this boundary condition is 
imposed on the propagator for the normal Asian option, using the method of images. Finally 
in Section HVl we draw conclusions. 



II. GEOMETRIC ASIAN OPTIONS IN THE BLACK-SCHOLES MODEL 



A. Partitioning the set of all paths 



The path integral propagator is used in financial science to track the probability distribu- 
tion of the logreturn x t = log (St /So) at time t, where So is the initial value of the underlying 
asset. This propagator is calculated as a weighted sum over all paths from the initial value 
Xq = at time t = to a final value xt = log(SV/So) a ^ time t = T : 

f-T 



K(x T ,T\0,0) = j Vx exp ^- C BS [x(t)} <ttj 
The weight of a path, in the Black-Scholes model, is determined by the Lagrangian 

on i2 



(1) 



Cbs W)\ 



1 

2^2 



x-\n 



(2) 
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where fi is the drift and a is the volatility appearing in the Wiener process for the logreturn 

For Asian options, the payoff is a function of the average value of the asset. Therefore 
we introduce xt = log (St /So) as the logreturn corresponding to the average asset price at 
maturity T. When St is the geometric average of the asset price, then xt is an algebraic 
average. 



x T 



T 



x(t)dt. 



(3) 



The key step to treat Asian options within the path integral framework is to partition the 
set of all paths into subsets of paths, where each path in a given subset has the same average 
xt- Summing over only these paths that have a given average xt defines the conditional 
propagator /C (xt, T |0, 0| xt)'- 

K(x T ,T \0,0\x t ) = Jvx6^x T -^^ x(t) dtj exp ^- £ BS {x(t)]d?J (4) 
This is indeed a partitioning of the sum over all paths: 



00 

/C(x T ,T|0,0) = J dx T JC(x T ,T |0,0|x T ) 



(5) 



The delta function in the sum J T>x over all paths picks out precisely all the paths that will 
have the same payoff for an Asian option. 

The calculation of /C (xt, T |0, 0| xt) is straightforward; when the delta function is rewrit- 
ten as an exponential, 



K(x T ,T 10, 01 x T ) 



dk c ikx 
2ir 



Vx exp (- J (c BS [x(t)) + ^ikx{t)^j dtj , (6) 



the resulting Lagrangian is that of a free particle in a constant force field in ID. The resulting 
integration over paths is found by standard procedures 24j : 



K(x T ,T |0,0|x r ) 



73 

7T(T 2 T 
6 



1 



exp 



2a 2 T 



x T — \ ^ 2 



xt 

a 2 T v XT " y 



(7) 



and corresponds to the result found by Kleinert [6j and by Linetsky 20]. 
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B. Link with stochastic calculus 



The conditional propagator K. (xt, T |0, 0| xt) is interpreted in the framework of stochas- 
tic calculus as the joint propagator K, (xt, Xt, T |0, 0, 0) of xt and its average Xt- The 
calculation of fC (xt, Xt, T |0, 0, 0) here is similar to the derivation presented in Ref. [25 1 
where this joint propagator is calculated for the Vasicek model. The main point is that in a 
Gaussian model the joint distribution of the couple {xt, xt} has to be Gaussian too. As a 
consequence this joint distribution is fully characterized by the expectation values and the 
variances of xt and xt and by the correlation between these two processes. The expectation 
value of Xt (t) is given by (ji — ^-j ~, its variance by ^ and the correlation between the 
two processes by The density function of such a Gaussian process is then known to be 



1C(x T ,x T ,T\0,0,0) 



x/3 
7ia 2 T 



exp 



X T 



a 2 T 

2 



X T 



a -r\T 



-i 2 



a 



Xt 



T 

L 1 ?r I TT 



Xt 



a 



T 
~2 



(8) 



This agrees with Eq. (j7]) for /C (xt, T |0, 0| xt)- 



C. Pricing of an average strike geometric Asian option 

If the payoff at time T of an Asian option is written as V^ smn (xr, Xt), then the expected 
payoff is 

oo oo 

E [V^^ixT, x T )] = J dx T J dx T V^ sian {x T , x T )K (x T , T |0, 0| x T ) (9) 

— oo — oo 

The price of the option, V^ sian is the discounted expected payoff, 

yAsian = e -rT E Asian ^ ^ 

where r is the discount (risk-free) interest rate. Using expression ((9]) the price of any option 
which is dependent on the average of the underlying asset during the lifetime of the option 
can be calculated. We will now derive the price of an average strike geometric Asian call 
option explicitly. In order to do this, expression (jUJ) has to be evaluated using the payoff: 

V^^ixT, x T ) = max(S T - S T , 0) = S max(e XT - e^ T , 0) (11) 
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Substituting (TTTTl in (TO yields 



V 



Asian 




S e 



-rT 



dx T / dx T (e XT - e* T ) K (x T , T |0, 0| x T ) 



(12) 



where the lower boundary of the xt integration now depends on xt- When considering an 
average price call, the payoff (for a call option) is max(SV — K, 0) leading to a constant 
lower boundary log (K/ So) for the xt integration, and the integrals are easily evaluated. 
In the present case however, the integration boundary is more complicated and it is more 
convenient to express this boundary through a Heaviside function, written in its integral 
representation: 

1 



V 



Asian 




oo r+oo gi(xy— iy)r 

dxx I dxx / dr : — 

oo J — oo 



T — IE 



(13) 



x (e XT -e* T )lC{x T ,T \0,0\x T ) 



Now the two original integrals have been reduced to Gaussians at the cost of inserting a 
complex term in the exponential. Expression ( f!3l) can be split into two terms denoted I\ 
and 1 2, where 



h = S e 




+oo 



dxi 



+oo 



dxi 



+ % {xt — xt) t + xt 



(14) 



a 2 T 2 J 

and I2 has the same form, except with xt instead of xt in the last term of the argument of 
the exponent. As a first step, the Gaussian integrals over xt and xt are calculated, yielding 

f(r) 



2?r« J_ ( 



-dr 



r — ie 



with 



/ (r) = exp 



a 2 T 2 / a 2 

r 2 + [a H 

6 V 2 



iT 



(15) 



(16) 



Now the integral has been reduced to a form which can be rewritten by making use of 
Plemelj's formulae. Taking into account symmetry, this reduces to 



fif) 



-dr = in 



T — IE 



erf 



2v^ 



1 



(17) 



with 



a 2 T 
6 



°1 \ Z 

Y) 2 



The first term thus becomes 

h = S e~ rT 



V - ^ erf 



na 2 T 2 



3T 
8^ 



/i + Y 



+ 1 



The second term, I2, is evaluated similarly, leading to 

Asian r» „—rT 

— ooe 








exp 



Using the cumulative distribution function of the normal distribution 



1 + erf 



x 

7! 



this can be rewritten in a more compact form as 



t rAsian c 

Vq — Joe 



with the following shorthand notations 

di = 
do = 



3T 




a 2 








3T 




a 2 


4a 2 


6 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



Expression (12"2"|) is the analytic pricing formula for an average strike geometric Asian call 
option, obtained in the present work with the path integral formalism. To the best of our 
knowledge, no pricing formula of this simplicity exists. To check this formula, we compared 
its results to those of a Monte Carlo simulation. The Monte Carlo scheme used is as follows 



251 ] : first, the evolution of the logreturn is simulated for a large number of paths. This 
evolution is governed by a discrete geometric Brownian motion for a number of time steps. 
Using the value for the logreturn at each time step, the average logreturn can be calculated 
for every path. Subsequently the payoff per path can be obtained, which is then used to 
calculate the option price by averaging over all payoffs per path en discounting back in time. 



The analytical result and the Monte Carlo simulation agree to within a relative error of 
0.3% when 500 000 samples and 100 time steps are used. This means that our analytical 
result lies within the error bars at every point. We also obtained the result for an average 
price Asian option; in contrast to the new result for the average strike option this could be 



compared to the existing formula 20|, |23( , and was found to be the same. 



III. ASIAN OPTION WITH A BARRIER ON A CONTROL PROCESS 
A. Derivation of the option price 

In this case we consider two stochastic processes: 

( ° 2 \ 
dx = l/i — — J at + adWi 

) J (24) 
dy= lu-^-Jdt + ^dZ 

which are correlated in the following manner: (dWidZ) = pdt. The x process models the 
logreturn of the asset price which underlies the Asian option, and the y process describes 
the control process. The payoff for an Asian option with a barrier on a control process 
is the same as for a normal Asian option, with the extra condition that the payoff is zero 
whenever the value of y surpasses a certain predetermined barrier. This is an example of an 
up-and-out barrier. There are other types of barrier options, namely down-and-out etc., but 
since their treatment is analogous we will not consider them here. The payoff for an Asian 
option with a barrier on a control process is given by: 

Vfi*T,*T) J S ^ T ~ K =»(«><»» (25) 

[ 3i<= [0,T] :y(t) >y B 

where the payoff of an average price Asian option has been used. Here Sq x denotes the initial 
asset price of the asset corresponding to the logreturn x and ys is the value of the barrier 
which has been placed upon the y process. It is difficult to price this option using payoff ( 125]) 
because of the extra barrier condition. However, if this condition could somehow be included 
in the propagator for these two processes, then the payoff would reduce to that of a normal 
(average price) Asian option, making the calculations more tractable. To construct this new 
propagator, henceforth called barrier-propagator ; a linear combination of propagators for the 
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combined evolution of both processes given in ([21]) can be taken: 

JCy (x T , y T , x T , T |0, 0, 0, 0) = K (x T , y T , x T , T |0, 0, 0, 0) + C K. (x T , y T , x T , T \x s , ys, x s , 0) 

(26) 

where fCi 3 stands for the propagator for the processes x and y where a barrier condition has 
been placed upon the y process. The propagator K (xt, yr, %t, T |0, 0, 0, 0) belonging to the 
system f l2~^|) is an extension of the propagator (J7j), and is given by: 



K (x T ,y T ,x T ,T |0, 0,0,0) 



x exp 



27r 3 T 3 a 4^2 (J _ p 2 
P 



1 



<7 



2a 2 (1 -p 2 )T 



<7 



25 2 (i-p 2 )r 



X 



2/T 



a 2 T 



2 - 




(27) 



Furthermore C is a factor upon which three conditions will be placed and xs,ys represent 
the initial condition from which the mirror-propagator starts. This mirror-propagator is 
used to eliminate all paths that cross the barrier, and because the paths represented by the 
mirror-propagator usually have higher values than the paths represented by the propagator 
/C (xt, Vt, xt, T |0, 0, 0, 0), they have been given another average x$- The barrier-propagator 
must be zero at the boundary: 



K* (x T ,y B ,x T ,T |0, 0, 0, 0) = 



(28) 



Using this boundary condition, an expression for C can be derived which must satisfy three 
conditions: firstly C must be independent of the averages xt and x$, secondly it may not 
depend on xt and finally it must be time- independent. This eventually leads to the following 
propagator for the total system of correlated stochastic processes x and y, with a barrier 
condition on y when y T G [—00, ys[- 



{xt, y T , x T , T |0, 0, 0, 0) = K (x T , y T , x T , T |0, 0, 0, 0) 



2 VB 
.€(4-3p 2 ) 



— 



JC(x s ,ys,x s ,T\0, 0,0,0) (29) 
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with the following shorthand notations: 

2y B pa 



t (4-3p 2 ) 
- (x s - X T ) + 



xs - x T ) 2 + 4 x T (x T - x T ) 



(30) 



k Vs = 2y B 



The propagator (1251) is equal to zero when y T G [y^,+oo]. A graphical presentation of 
propagator (1291) is shown in Fig. [TJ 




Figure 1: Graphical representation of the barrier-propagator (|29p in arbitrary units, for the system 
of two correlated processes x and y, given by (j24[) . where a barrier has been placed on the y 
process at ys = 0.2. The following values were used in this figure: = v = 0.03 Year- 1 ;^ = £ = 



0.25 Year _ 2;T = 1 Year 
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Using the propagator (I29I) the price of an Asian option with a barrier V^ B can be calcu- 
lated. The general pricing formula is given by: 



V AB = e~ rT 



dxi 



dx T K B (x T , yr, x T , T |0, 0, 0, 0) V™ an (x T , x T ) (31) 



This calculation was done for an average price option: Vf Lsian (xT, xt) = max (So x e XT — K, 0). 
The calculation, though rather cumbersome, is essentially the same as for the Asian options 
in section [III The integral over xt is a Gaussian integral, and the remaining two integrals 
can be transformed into a standard bivariate cumulative normal distribution, defined by: 

«6 



N[a, b; x\ 



1 



27r v /r^ 



exp 



oo J — oo 



2(1-X 2 



(x 2 + y 2 ~2 X xy) dxdy (32) 



This eventually leads to the following pricing formula for an Asian option with a barrier: 



V, 



AB 



-rT 



N I d 1 ,d 2 , -\/ -p 



K N I d 3 ,d 4 , 



3T ( X S 



■)[■*+(<-#)] (*- 

S y 



«(4-3p 2 ) 



iV I d 5 ,d 6 , 



V Si 



Oy 



e(4-3 P ^) 



/ 



iV I d 7 ,d 8 , -p 



(33) 
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where the following shorthand notations were used: 



In 



di 



K 

Sqx 



T ( a 2 
2 ^1 



In 



d 2 



B 

S y 



a 2 T 



e 

2 



1 SoxJ 2 V 2 



In 



B 



a 



2 T 



In 



K 

Sqx 



T 



2 T + h " 



a 
~6 



2\ 1 



In 



da 





T 


(soy) 


a 



3£p- + a- + a^- y + — 



T 



2xs 1 f o 2 



In 



d* 



B 



O.y 



T 
27 



c T 



- (6p^ + 2ay s ) + 2a(u-^ 



(34) 



B. Results and discussion 



Fig. ([2]) shows the option price for an Asian option with a barrier as a function of the 
initial asset price belonging to the y process, defined by: Sr y = So y e yT . This figure shows 
that the analytical result derived in section II II Al deviates from the Monte Carlo simulation 
with increasing correlation. The approximate nature of our approach can be understood as 
follows. The essence of the approach presented here is that to calculate the price of Asian 
barrier options, two steps need to be taken. First, a partitioning of paths according to the 
average along the path must be performed, and second, the method of images must be used 
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Figure 2: Option price for an Asian option with a barrier, as a function of the initial asset price 
belonging to the y process. Since no specific asset is considered the option price could be stated in 
any currency, therefore this figure is given in arbitrary units. The analytical result deviates from 
the Monte Carlo simulation (MC) for increasing correlation. The value for the barrier used in this 
figure is B = 150. 



in order to cancel out paths which have reached the barrier. The difficulty combining these 
two steps, is that mirror paths have a different average than the original paths, and thus 
belong to a different partition. This difficulty can apparently be overcome by treating the 
average itself as a separate, correlated process (as proposed in Ref. Q). This procedure, 
relating /C (xt, T |0, 0| xt) to K (xt, %t, T |0, 0, 0), leads to the correct propagator (and price) 
in the case of a plain Asian option as shown in section III Bl 

However, from the results shown in Fig. [2] it is clear that this is no longer the case for an 
Asian option with a barrier on a correlated control process. This is because the exact average 
of the x process does not behave as a separate, correlated process (the average described 
by this process is henceforth called the approximate average). This approach is exact for a 
plain Asian option, where all paths contribute, but when a barrier is implemented using the 
method of images, and thus eliminating some of the paths, the following approximation is 
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made. When the y process hits the barrier and is thus eliminated, its corresponding x and 
x processes are eliminated as well. But the x process considered in our derivation is only 
approximate, so the wrong x paths are eliminated. The central question is whether this will 
lead to a difference between the distribution of contributing paths for the exact averages 
and the corresponding distribution for the approximate averages, when a barrier has been 
implemented. Figure 3 shows that this is indeed the case, and that this difference increases 
when correlation increases. When the correlation is zero, the paths which are eliminated for 




Figure 3: Comparison of the distribution of the exact and the approximate average of the logre- 
turn x, when a barrier is present, a) When the correlation equals zero, both distributions are equal 
b) For relative small values of the correlation, the assumption that the x process is correlated with 
the process for the logreturn x is still relatively good, c) When the correlation is relatively large, 
the distribution of the exact averages deviates significantly from the approximate averages. 

both the exact and the approximate average are randomly distributed (because the behavior 
of x has nothing to do with the behavior of y) , which means that both distributions remain 
the same Gaussian as they would be without a barrier. This is the reason why our result 
is exact when correlation is zero. Another source of app roximation lies in the use of the 
Black-Scholes model which has well-known limitations [lj, [l^] ■ Several other types of market 
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models propose to overcome such limitations, for example by introducing additional ad hoc 
stochastic variables {26] or by improving the description of the behavior of buyers/sellers 
[27)]. The extension of the present work to for example the Heston model lies beyond the 
scope of this article. 



IV. CONCLUSIONS 



In this paper, we derived a closed-form pricing formula for an average price as well as an 
average strike geometric Asian option within the path integral framework. The result for the 
average price Asian option corresponds to that found by Linetsky 20], using the effective 



classical partition function technique developed by Feynman and Kleinert 15]. The result 
for the average strike Asian option was compared to a Monte Carlo simulation. We found 
that the agreement between the numerical simulation and the analytical result for an average 
strike Asian option is such that they coincide to within a relative error of less than 0.3 % 
for at least 500 000 samples and 100 time steps. 

Furthermore, a pricing formula for an Asian option with a barrier on a control process 
was developed. This is an Asian option with the additional condition that the payoff is 
zero whenever the value of the control process crosses a certain predetermined barrier. The 
pricing of this option was performed by constructing a new propagator which consisted of 
a linear combination of two propagators for a regular Asian option. The resulting pricing 
formula is exact when the correlation is zero, and is approximate when the correlation 
increases. The central approximation made in our derivation, is that the process for the 
average logreturn x is treated as a stochastic process, which is correlated with the process of 
the logreturn x. This assumption is correct whenever all price-paths contribute to the total 
sum, but becomes approximate when a boundary condition is applied. 
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